arXiv:1505.05282vl [hep-ph] 20 May 2015 
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We perform a dispersive analysis of the ljtt electromagnetic transition form factor, using as input 
the discontinuity provided by unitarity below the lutv threshold and including for the first time 
experimental data on the modulus measured from e'*"e“ —>■ wtt® at higher energies. The input 
leads to stringent parameterization-free constraints on the modulus of the form factor below the 
LOTT threshold, which are in disagreement with some experimental values measured from uj —>■ 7r°7* 
decay. We discuss the dependence on the input parameters in the unitarity relation, using for 
illustration an N/D formalism for the P partial wave of the scattering process lott ^ tttt, improved 
by a simple prescription which simulates the rescattering in the crossed channels. Our results confirm 
the existence of a conflict between experimental data and theoretical calculations of the um form 
factor in the region around 0.6 GeV and bring further arguments in support of renewed experimental 
efforts to measure more precisely the ui decay. 


I. INTRODUCTION 

The transition form factors of light mesons play an im¬ 
portant role in low energy precision tests of QCD [I| . In 
particular, they enter as contributions to hadronic light- 
by-light scattering calculations 0 , which are crucial for a 
more accurate theoretical determination of the standard 
model prediction for the muon’s anomalous magnetic mo¬ 
ment (for recent reviews see i, i). 

The case of the ojtt electromagnetic form factor is par¬ 
ticularly interesting as there are some discrepancies be¬ 
tween the theoretical calculations and the experimental 
data from the decay w —>■ 7r°7* reported in @-13 • This 
form factor was described by Vector Meson Dominance 
rVMD) model and by a chiral Lagrangian approach in 
[ 1 , @1 . Calculations based on a standard dispersion re¬ 
lation were performed a long time ago in [ij and re¬ 
cently in . The discontinuity of the form factor 

required in the Cauchy integral can be expressed in terms 
of known observables by using unitarity. The two-pion 
contribution to the unitarity sum gives the discontinu¬ 
ity in terms of the P partial wave of the amplitude of 
the process ojtt — > tttt, itself calculated in the dispersion 
theory Emm, and the pion electromagnetic form 
factor, a quantity which is known with very good preci¬ 
sion. However, the two-pion approximation is valid only 
in a region which extends to a good approximation up 
to the LOTT threshold, = (m^ -|- m^)^. Due to the lack 
of information on the discontinuity above this threshold, 
various assumptions were adopted for the evaluation of 
the dispersion integral, either by appl ying th e two-pion 
approximation also at higher energies [2, [ 13 , or by ex¬ 
panding the dispersion integral in powers of a suitable 
variable [l^ . 

A study performed recently in [l^ used as input above 
the OJTT threshold, instead of the discontinuity, a model- 
independent integral condition on the modulus squared 
of the OJTT form factor. The condition was obtaine d by 
using an approach proposed originally by Okubo [iSj, 
which has come to be known as the method of unitar¬ 


ity bounds (a recent review of this approach is presented 
in [3). It exploits unitarity and the positivity of the 
spectral function of a suitable current-current correla¬ 
tor, calculated by operator product expansion (OPE) in 
the Euclidean region. In the particular case of the cjtt 
form factor, the method, adapted to the specific input 
conditions available, led eventually to a functional opti¬ 
mization problem of a type considered for the first time 
in [Ulii- The solution of the problem yields upper and 
lower bounds on the modulus of the wtt form factor in the 
region 4m^ < t [3|. A specific feature of this form 

factor is that its discontinuity across the cut is not purely 
imaginary. As a consequence, the form factor is not a 
real analytic function, as happens in familiar cases like 
the pion vector form factor. Therefore, in M the formal¬ 
ism of bounds was extended to functions which are not 
real analytic. Although not very stringent, the bounds 
derived in [l3| are in disagreement with the experimen¬ 
tal data on the modulus of the form factor in the region 
around 0.6 GeV, measured from the decay w —7r°7*, 
confirming thus the conclusion of the analysis [ll| based 
on standard dispersion theory. 

It is important to note that the dispersive analyses per¬ 
formed so far did not include experimental data on the 
form factor available in the scattering region, above the 
WTr threshold. Measurements of the modulus from the 
reaction e+e“ —>■ are reported in [Tol - I^ (a set of 

such data is shown in Fig. [TJ where we show for com¬ 
pleteness also the modulus measured in the decay region 
t < t-, t- = {rjiuj — m,n.)^). In the present paper we 
consider the problem of including this information in the 
dispersive formalism. Specifically, we perform an anal¬ 
ysis of the form factor using as input the discontinuity 
for t < t+, calculated in a theoretical model based on 
unitarity, and experimental information on the modulus 
for t > Even though the modulus is not known at 
all energies, we can implement the information in a con¬ 
servative way, as a condition on a weighted integral of 
the modulus squared from to infinity. This leads to 
a mathematical problem similar to that encountered and 
solved in M- The result is expressed in the form of ex- 


plicit upper and lower bounds on the modulus of the wtt 
form factor below calculable in terms of the disconti¬ 
nuity below 4 and the modulus above . The formalism 
provides therefore a consistency test for the experimen¬ 
tal data on the wtt electromagnetic form factor, which 
exploits analyticity and unitarity in a parametrization- 
free way. 

The theoretical input of the test consists from the uni¬ 
tarity relation giving the discontinuity of the form factor, 
which involves the amplitude of the ttw —>■ tttt scattering. 
It is of interest to study the influence of the possible 
uncertainties in this part of the input. The full calcu¬ 
lation involves the solution of integral equations known 
as Khuri-Treiman equation s 12611. The amplitudes ob¬ 
tained in this formalism in |12l Il3l| are available only in 
numerical form. The older treatment performed in [lOj, 
based on N/D formalism, has the advantage of providing 
explicit expressions. However, it is not entirely satisfac¬ 
tory, as it does not account for the rescatterings between 
all the final pions in the kinematical region where the 
decay w —>■ tttttt is allowed. It is worthwhile trying to 
cure the shortcomings of this approach. In this paper 
we consider an improved N/D treatment, obtained by 
applying to a prescription proposed in for in¬ 
cluding finite-width effects in the resonances exchanged 
in the crossed channels. The improved model captures 
the essential features of the full solution, still preserving 
the explicit dependence on the input parameters. This 
has enabled us to investigate the influence of the uncer¬ 
tainties of the theoretical part of the consistency test. 

The paper is organized as follows. In the next section 
we review the basic definitions and show how the input 
information on the form factor can be expressed as an 
extremal problem for analytic functions. In section IIIII 
we review the main steps of the proof and write down 
the solution of the extremal problem obtained in Ref. 
H- The results are presented in Sec. El where we 
discuss also their dependence on the parameters of the 
input. Sec. El contains our conclusions. The paper has 
an Appendix where we describe briefly the N/D model 
of [1^ and an improved version based on a prescription 
suggested in 


II. INPUT IN THE CONSISTENCY TEST 


We begin with a brief description of the form factor and 
the constraints that it satisfies. We use the conventions 
of where the form factor funr{t) is defined from the 
matrix element 


(w(pa, A)7r°(p6)|j^(0)|0) = 

( 1 ) 

where 4 is the isovector part of the electromagnetic cur¬ 
rent, A denotes the w polarization, q = Pa+Pb and t = q^. 
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FIG. 1: Modulus of the oj-k form factor measured from 
ui 7r°7* decay by Lepton-G @], NA60 (2009) and 
NA60 (2011) [3, and from e^e“ ^ oitt® by GMD-2 (2005) 

M- 


In the convention adopted here^ the form factor funrit) 
has dimension of GeV“^. 

Unitarity implies that fui-wit) has a cut along the real 
axis for t > 4m^. Keeping the two-pion contribution in 
the unitarity sum, the discontinuity of fujTr{t) across the 
cut is written as 

discfajnit) = ^^F*(t)fi{t)9{t-Aml), t<t+, (2) 


where q{t) = \Jt/d — is the center of mass momen¬ 
tum of the pion pair, is the pion electromagnetic 

form factor and fi (t) the P partial-wave amplitude of the 
scattering 


w(pa, A)7r°(pb)-)> 7r+(gi)7r (ga)- (3) 

The scattering process is physical for t > t+. In the 
region 4m^ < t < t-, where the decay uj —>■ 7r+7r“7r° 
is allowed, /i(t) is the P-wave projection of the decay 
amplitude, while the region < t < t_|_ is unphysical. 

The amplitude fi {t) was calculated in in the frame 
oi N/D formalism, with the left-hand cut described by 
poles in the crossed channels due to the exchange of the 
p meson assumed to be stable. In this model, the phase 
of fi(t) coincides with the tttt P-wave phase shift and 
exactly compensates in the discontinuity ([2) the phase of 
F*{t), related also to the tttt P wave phase shift by 


^ The dimensionless form factor FT^u;-y(t) defined in IlCll is re¬ 
lated to the definition adopted here by A-n-(t). 

The form factor F-^uiit) defined in 1271 is dimensionless, nor¬ 
malized as FTrui(O) = 1 and is related to our definition by 
fuin(t) = 2 \/2Cb Uru; (t), where Cu} is defined in Eq. (35) of 
\21V in terms of the total width 
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Watson theorem [2^. Therefore, the form factor funr{t) 
calculated in is a real analytic function^. 

In the more complete calculation [IIHIl based on 
Khuri-Treiman formalism [2^ , the amplitude /i (t) is ob¬ 
tained by numerically solving a set of integral equations. 
Due to the rescattering between the final pions in the 
decay region 4m^ < t the phase of fi{t) does not 

coincide with the tttt P-wave phase shift, as one would 
naively expect from Watson theorem. Therefore, the 
phases of the two factors in (ED do not compensate each 
other, and the discontinuity is not purely imaginary 
[TJ, [l3|. In other words, the ojtt form factor is not a real 
analytic function, which is true also in the case of other 
transition form factors (^ . 

The expression @ is valid only in the region 4m^ < 
t < t+, since above the con threshold other intermedi¬ 
ate states, besides the two-pion states, contribute in the 
unitarity sum. So, strictly speaking the discontinuity of 
fuiTT (t) is not available for t > On the other hand, the 
modulus of fojTvit) can be extracted from experimental 
data on e+e“ —process, measured in The 

connection between the cross section and the modulus 
[Sill is, in our convention, 

= (4) 

where p(t) = {t — t-){t — t+)/At is the center of mass 
momentum of the con pair in the rest system of the virtual 
photon and we recall that t± = (m^ ± 

Using the experimental data on the modulus and the 
asymptotic behaviour \f^{t)\ ~ 1/t predicted by per¬ 
turbative QCD scaling [^, it is possible to obtain a rea¬ 
sonable estimate of a weighted integral over the modulus 
squared from t+ to infinity. Thus, we can write an L^- 
norm condition of the form 

1 r°° 

- |^7r(t)P = /, (5) 

Jt+ 

where w{t) is a suitable weight, chosen such as to allow 
a precise evaluation of the quantity I. We mention that 
a similar way of including experimental information on 
the modulus at higher energies was adopted in recent 
investigations of the pion electromagnetic form 

factor. 

In the present analysis we have considered weights of 
the simple form 


w{t) = 


1 

U’ 


( 6 ) 


^ A function F{t) analytic in the t-plane cut for t > 4m^ is of real 
type if it satisfies the condition F{t*) = {F{t))*. In particular, 
this implies that the function is real on the real axis for t < 4m^, 
and its discontinuity across the cut can be written as disc F(t) = 
F{t + ie) — F(t — it) = 2i Im F{t + it). 


where the value of c > 0 is taken such as to suppress the 
high-energy tail of the integral, where the form factor is 
not known. In practice we evaluated the quantity / using 
an interpolation of the data on modulus from [22| shown 
in Fig. [1] from t+ = 0.84 GeV^ up to t = 1.86 GeV^, 
continued in a smooth way with a modulus \fu-Kit)\ de¬ 
creasing like 1/t. As will be clear in the next section, for 
a fixed weight the results of the formalism depend mono- 
tonically on the numerical value of /, in the sense that a 
larger I gives weaker results. Therefore, for a conserva¬ 
tive estimate, we have used as input in the data region 
the central values from enlarged by their quoted er¬ 
rors. For c = 2 this leads to 

I = 4.63 GeV“^ (7) 

where the region above 1.86 GeV^ contributes to the in¬ 
tegral with about 8 %. 

Finally, we use as input the value of |/t.j7r(0)|, known 
experimentally from the uj —> decay rate. The up¬ 

dated value is [3^ 

|/^,(0)| = (2.30 ±0.04) GeV-^ (8) 

As already mentioned, in the present paper we shall 
check the consistency of the data shown in Fig. [T] by 
comparing the data in the decay region with the allowed 
range of the modulus imposed by unitarity and analyt- 
icity. Mathematically, the problem amounts to deriving 
upper and lower bounds on \fanT(t)\ for t < t+, upon 
the class of functions funv(t) analytic in the t-plane cut 
for t > 4mwhich satisfy the following conditions: {i ) 
their discontinuity is given by in the region t < t+^ 
(ii) they satisfy the constraint and {Hi) they sat¬ 
isfy the condition ([5]). The solution of this mathematical 
problem will be given in the next section. 


III. SOLUTION OF THE EXTREMAL 
PROBLEM 

An extremal problem of t he ty pe mentioned above was 
solved for the first time in [13, [lal on the class of real 
analytic functions. The generalization to functions which 
are not real analytic was investigated in detail in M- We 
do not repeat the whole proof here, but only outline the 
main steps and write down the solution. 

The first step is to map the t plane cut along t > 4 
onto the unit disk \z\ < 1 in the z = z{t) plane. We have 
adopted the conformal mapping 


z{t) 


1 - yi - t/t+ 

1 ± \/l — t/t+ 


(9) 


which brings the origin of the t plane to the origin of 
the z plane, z(0) = 0. In the z-plane the elastic region 
4m^ < t < becomes the segment < x < 1 of 
the real axis, where x^n- = 5(4m^) > 0, and the upper 
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(lower) edges of the cut t > t+ become the upper (lower) 
semicircles^. 

Further, we construct a so-called outer function [s^, 
i.e. a function analytic and without zeros in \z\ < 1, its 
modulus on \z\ = 1 being equal to \/w{i{z))\di{z)/dz\, 
where w{t) is the weight appearing in ([5]) and i{z) is the 
inverse of the function z{t) defined in ([5|). The general 
expression of the outer functions is given in (see also 
the review 0). For weights w{t) of the form dHl), we 
obtain for the outer function, denoted as C{z), the exact 
analytic expression [T^ 

C{z) = (2v^)i-"(l - z)i/2(i + (10) 

From this expression it follows that C{x) is real and pos¬ 
itive for — 1 < a: < 1, which corresponds in the t-plane to 
the semiaxis t 

If we introduce now a new function h(z) by 

h{z) = C{z) f^^{i{z)), (11) 

the condition ([5|) takes the simple form 

d0\h{e^^)f = I. (12) 

The function h(z) is analytic in \z\ < 1 except for a cut 
along the segment (x^, 1), where its discontinuity is 

disc h{x) = A{x) = C{x) disc funr{i{x)). (13) 
By expressing h{z) as 

h{z) = ^[ ^^dx + g{z), (14) 

27rj X — z 

the new function g{z) is analytic in \z\ < 1, as its discon¬ 
tinuity across the cut vanishes: 

disc (?(a;) = 0, a:€(a;^,l). (15) 



Since we consider in general form factors that are not 
real analytic, the function g(z) is analytic, but its values 
on the real axis may be complex. 

We now express the available information on the form 
factor as a number of constraints on the function g. By 
inserting (HI in dH we obtain the condition 


1 

27r 




A(x) 

r» _ 


dx + (?(e*®) 



(16) 


and using m and m we write g{0) as^ 




^.(O)C(O) 


1 

2TTi 


A(x) 


dx. 


(17) 


^ Other mappings are obtained by changing the point that is 
mapped onto the origin of the 2 :-plane. It can be shown [T3l 
that the results do not depend on the choice of the conformal 
mapping. 

^ Note that there is a misprint in the expression of ^(0) given in 
Eq. (24) of Ref. jjj. It did not affect the results of [l^ since 
the calculations were peformed with the correct expression. 


The problem is to find the maximal allowed range of 
|^(zi)| at an arbitrary given point zi = z{ti) in the in¬ 
terval for functions g{z) analytic in < 1 and 

subject both to the boundary condition (fT5|) and the ad¬ 
ditional constraint (HZD- 
It is useful to denote 


9{zi) = C, (18) 

where ^ is an unknown parameter. Then one can prove 
(see for instance [l^) that the allowed range of ^ is de¬ 
scribed by the inequality® 

m 2(C < r (19) 

where is the solution of the functional minimization 
problem 


dliO 


mm — 
2tt 



1 

2TTi 


A{x) 


dx + g(e*®) 


X — e 


id 


( 20 ) 


upon the class of functions analytic in | 2 :| < 1, 
which satisfy the constraint (Ha and the additional con¬ 
dition (IT51) for a given 

The constrained minimum norm problem (EHl) was 
solved in [l^ by the technique of Lagrange multipliers, 
leading to a solution written in compact form: 


dl{i) 


1 

47r^ 


A{x)A*{y) 

l-xy 


dxdy -h |5(0)l^ 


Tie-<?(0)r (21) 

Zi 


By inserting (EH) in EH) we obtain upper and lower 
bounds on the parameter Expressed in terms of the 
form factor funr{t) by using Eqs. (fTTl) and (fUl) . they lead 
to the inequalities |1^ : 


\fu,At)\ < 


ff(9) + ^ Ix 


x — z{t) 


dx 


3(0) 


1 

2771 JX 


cm) 

1 A(3;) 


x — z{t) 


dx 


cm) 




id)i' 


( 22 ) 


where z{t) G {xt^,1) is the image of the point t in the 
z-plane, and 



A(x)A*(^) 

l-xy 


1/2 


dxdy - |g(0)p 


(23) 

We recall that C{z{t)) > 0 for t < t+, which justifies its 
appearance outside the modulus sign in the denominator 
of (1221). 


^ This shows that the results remain the same if 0 is replaced by 
an inequality involving a quantity that majorizes I. 






























The upper and lower bounds (E21) are calculable in 
terms of the input defined in the previous section. They 
determine an allowed interval for the modulus \funr(t)\ at 
every t < t+. From (1^^ it follows that, for a fixed weigth 
w{t) in the T^-norm constraint ([5]), the bounds depend 
monotonically on the value of I : smaller values of / lead 
to narrower allowed intervals for \funr{t)\ at t < t+. We 
already took into account this property for a conservative 
estimate of I, as discussed in the previous section. 

It is useful to remark also that, since the last term in 
m is positive, from m and m we can write down 
the inequality 


1 

47r^ 


A{x)A*{y) 

1-xy 


dxdy+ |g(0)p < I, 


(24) 


where (?(0) is defined in (IT7)) and A(a;) is (fT^ . The 
inequality (HI involves only input quantities and rep¬ 
resents a necessary condition that must be satisfied by 
them. If it is violated, the input is not consistent with 
analyticity and unitarity. 


IV. RESULTS 

We have investigated several suitable weights of the 
form m and checked that they lead to similar results. 
The calculations reported below were done with the 
choice c = 2, which ensures a good suppression of the 
high energy part of the integral. 

As already mentioned, we have employed the discon¬ 
tinuity (H|) of the WTT form factor in the range (4rni,t-t-) 
from the recent dispersive treatment reported in [ll| and 
from the older work [l^. In m the pion vector form 
factor FV(i) has been reconstructed from an Omnes rep¬ 
resentation using as input the pion~pion phase shift 

(t) calculated from Roy equations in [H, . In 0, 

the pion form factor was described by a Gounaris-Sakurai 
representation given in Eqs. (lA. 1411 - (lA. 161) of the Ap¬ 
pendix. We have checked that the differences between 
the two representations of the pion form factor are very 
small and have a negligible influence on the results. 

On the other hand, the differences in the partial wave 
fi{t) used in M and m are sizable, and have a larger 
impact on the results. In Fig. [5] we show the allowed 
bands determined by the upper and lower bounds on the 
modulus squared (normalized to its value at t = 0) in 
the part of the elastic region accessible experimentally 
in a; —>■ decay, calculated using the expressions 

(1^^ with input from [1^ and [HI- The results shown in 
Fig. [2] were obtained by varying the input value at t = 0 
inside the error bar given in ([5]) and taking the weakest 
bounds, i.e. the largest allowed bands at each energy. 
For comparison, we also show the result of the dispersive 
calculation performed in [H, and several experimental 
data from [H-13- 

With the input discontinuity from [l^, the allowed 
band is consistent with the dispersion relation calculation 



FIG. 2: Upper and lower bounds compared with experimental 
data on \fuin{t)/fujn{0)f- Cyan band: bounds calculated us¬ 
ing in the discontinuity m the partial wave amplitude fi{t) 
from 0- Grey band: bounds calculated using in the dis¬ 
continuity the amplitude fi{t) from m- The yellow 
band is the result of the dispersive calculation performed in 
[lll| . The data are from Lepton-G 0], NA60 (2009) @| and 
NA60 (2011) [3]. 


performed in that work. The allowed band obtained with 
the partial wave fi{t) from [ifij i is shifted upwards and 
the two bands do not overlap. For both inputs the upper 
bounds shown in Fig. H] are significantly lower than the 
data from in the region around 0.6 GeV. 

We mention that the upper and lower bounds shown in 
Fig. [2] are much more stringent than the upper and lower 
bounds obtained in 0 with the same input on the dis¬ 
continuity m, but with a model-independent condition 
on the modulus above the ujtt threshold®. 

It is of interest to understand the origin of the differ¬ 
ence between the two predictions shown in Fig. [21 To 
this end, we have calculated the bounds using also the 
improved version of the N/D model of [l^, which in¬ 
cludes the effect of rescattering in the crossed channels as 
discussed in the Appendix. The improvement has the ef¬ 
fect of shifting the bounds downwards, towards the band 
calculated with the input amplitude fi (t) from [Tl|, [I^ , 
but the shift is small, of a few percents. For illustra¬ 
tion we present in Fig. [3] the bounds calculated with the 
improved N/D model for fi{t) in the whole range t < t+. 

As follows from Eqs. (IA.5I1 and (IA.19|) of the Ap¬ 
pendix, the N/D model contains as input the dimen¬ 
sionless coupling constants ypT^-rr and The results 


For instance, the allowed range of the ratio l.f^irfU/.A-n-fOlP at 
0.64 GeV obtained in H with input discontinuity from flUl is 
(1.6, 36.8), while the range predicted in this work with the same 
discontinuity is (11.3, 15.7). 
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FIG. 3: Orange band; bounds on \fujTrit)/fuj-n-(0)f in the 
whole region t < t+, obtained with the improved 
model. The data are from Lepton-G [^, NA60 (2009) [^, 
NA60 (2011) 0] and GMD-2 (2005) [13]. 



FIG. 4: Upper and lower bounds on |/ij7r(t)//u;7r(0)|^ calcu¬ 
lated using the improved N/D model, with gu}pTr = 11.5 (green 
band) g^pTr = 13.42 (orange band) and Quip-n = 15.5 (indigo 
line). The dashed lines are the bounds calculated with the 
discontinuity ® from 0. The red line is the form factor 
calculated in 0 [3| with a low-energy effective theory. 


presented in Fig s. [2] and [3] have been obtained using the 
values from 


5p777r = 5.96, = 13.42. (25) 

The quantity gp-j^Tr has not changed significantly over the 
last 40 years, the value (l25ll being fully consistent with 
the PDG-2014 width Fp, quoted below Eq. (IA.16|) . On 
the other hand, the coupling guipTv is still not very well 
known. We note that the dimensionless constant guipn 
used here is related to the similar parameter used in 
[ 23 , 1131 , which we denote as gippn to avoid confusion, by 


gcjpTT — Q uin-K - The values gupp-K — (13.8 i 0.3) GeV 
obtained in from a global fit of the wtt form factor, 
and 5 (pp 7 r = (16.7 ±0.4 ±0.6) GeV“^ derived in [23| from 
a fit of the cross section of e+e“ —>■ W7r°, correspond 
in our notation to the values gajp-K = 10.80 ± 0.23 and 
gujp-K = 13.07 ± 0.31 ± 0.47, respectively. We note that, 
while the above references consider parametrizations of 
the UJTT form factor, in the present analysis the coupling 
<?ipp 7 r is an input parameter in the model of the amplitude 

One might ask whether a suitable choice of the pa¬ 
rameter gpjpTT can reduce the conflict between the calcu¬ 
lated bounds and the experimental data. We have inves¬ 
tigated this question using the improved N/D model dis¬ 
cussed in the Appendix. We remark first that for values 
guip-TT > 15.5 the inequality (l24)) . which expresses a con¬ 
sistency condition on the input, is violated. Therefore, 
values of g^jp^ larger than 15.5 are not allowed, being in¬ 
consistent with the other input quantites and the general 
properties of analyticity and unitarity adopted here. 

In Fig. |3]we show the allowed bands calculated with 
three choices of the couplings guip-jr- the value (1331) . a 
smaller value, gujpTr = 11.5, and the maximum allowed 
value gpjpT^ = 15.5 mentioned above. In the latter case, 
when the the inequality (1331) is saturated, the quantity /' 
defined in (1231) is zero and the upper and lower bounds 
written in dUD become equal. Therefore, as shown in 
Fig. a the allowed band of the modulus shrinks in this 
case to a line. The curves show that the bounds ex¬ 
hibit a monotonous dependence on the value of gcpp^ the 
band obtained with g^jpp^ = 11.5 being consistent with 
the allowed domain obtained using as input fi(t) from 
the calculation [Hill- An important conclusion is that 
the disagreement with the experimental data around 0.6 
GeV is preserved even if the coupling g^pT^ is increased 
up to its maximum allowed value. 

We show also in Fig. [4] the ujtt form factor calculated 
in [3, I3| within an effective field approach. It exhibits 
a more rapid increase and slightly exceeds the highest 
allowed band shown in Fig. [4] near 0.6 GeV. 


V. DISCUSSION AND CONCLUSIONS 

Our study has been motivated by the existence of 
certain discrepancies between the recent calculations 
0113 , [11 of the LOTT electromagnetic transition form 
factor in the frame of dispersion theory, and the data 
measured from the uj —>■ 7r°7* decay around 0.6 GeV. 
The present work differs from the previous calculations 
by the different input used above the ujtt threshold 
while the investigations [ll|^ HI are based on a standard 
dispersion relation requiring the knowledge of the dis¬ 
continuity of the form factor along the whole cut, and 
the work 0 exploits a model-independent integral con¬ 
dition on the modulus, derived from unitarity and per¬ 
turbative QGD, we have resorted to experimental data 
obtained from e“'"e“ —> In order to reduce the bias 
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due to the absence of data at higher energies, we have 
implemented this information in a conservative way, as a 
weighted integral @ of the modulus squared. 

The aim of our study was to test the consistency of the 
experimental and theoretical information available on the 
WTr form factor in a parametrization-free approach. We 
have derived upper and lower bounds on the modulus for 
t <t+, using as input the discontinuity ([2|) in its region 
of validity below the wtt threshold, and the condition ([^ 
on the modulus above the um threshold. Mathematically, 
the problem is of the type considered some time ago for 
real-analytic functions in [Sill and generalized rece ntly 
to analytic functions which are not of real type in |14| . 
Since we used experimental data above the results ob¬ 
tained in the present work are much stronger than those 
obtained in (l4| . where only a theoretical inequality on 
the modulus above was exploited. 

The main theoretical ingredient of the analysis is the 
partial wave /i(t) entering the discontinuity (12). There¬ 
fore, it is of interest to establish the influence of various 
parameters entering this quantity on the final results. 
The amplitude fiit) calculated in [T^, E^l is available 
only in numerical form and is not suitable for this pur¬ 
pose. We have considered therefore the older calculation 
based on N/D formalism performed in [13 , which has 
the advantage of displaying in an explicit way the depen¬ 
dence on various parameters, and have improved it by a 
prescription suggested in for including the effect of 
rescattering in the crossed channels. 

Our study has showed that including the rescattering 
has the effect of shifting down the allowed band for the 
modulus of the form factor in the region t < . However, 

in the frame of the N/D model the effect is quite mod¬ 
est, of a few percents. On the other hand, the bounds 
are quite sensitive to the coupling g^ip-w, which enters as 
input in the calculation of fi{t) in the N/D formalism. 
The results obtained with /i(t) from can be re¬ 

produced by using the value guip-K = 11.5 in the improved 
N/D formalism. It turns out that values of g^ip-K larger 
than 15.5 are excluded, being inconsistent with the re¬ 
maining elements of the input. By increasing g^ipTr, the 
allowed bands are pushed upwards. However, as shown 
in Fig. 12 the narrow band calculated with the maximum 
allowed value of guipTr is still significantly lower than the 
experimental data from [i-Q near 0.6 GeV. 

Our results reveal a clear conflict between the experi¬ 
mental data on the modulus of the lott form factor mea¬ 
sured in the decay region t < t- from w —>■ 7r°7* and 
in the scattering region t > t+ from e+e —?► W7r°. We 
note that possible discrepancies between the data on the 
modulus measured at energies below and above 4 
can be noticed also in the attempts to describe the form 
factor with specific parametrizations (2^. In contrast, 
no parametrization of the form factor was necessary in 
the present analysis. The present work confirms the con¬ 
clusions of other recent dispersive analyses iHlil and 
brings further arguments in support of renewed experi¬ 
mental efforts to measure more precisely the uj conversion 


decays [3^. 1^. 
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Appendix: Improved N/D treatment of the ujn tttt 
amplitude 


The N/D model proposed in does not include the 
rescattering between all the final pions in the kinematical 
region where the oj decay to three pions is allowed. In 
this Appendix we briefly describe the model and present 
a simple modification, which is able to capture the char¬ 
acteristic features of the full solution. 

For convenience, we use in this Appendix the notation 
of M- The relation with the conventions used in the 
text is clear by comparing Eqs. (5.1) and (5.3) of M 
with Eqs. m and ((4]) of this paper, respectively. The 
P partial wave amplitude of the scattering process m, 
denoted in M as has dimensions of GeV ^ and is 
related to the partial wave /i (t) of pm by 

ti(t) = |m^/i(t). (A.l) 

In the N/D formalism, the amplitude t^{t) is written as 

m 

t^{t) =tl{t) (A.2) 

where t/{t) has only a left-hand cut and t/j(t) has only 
a right hand cut for t > 4m^. 

The piece t/{t) was calculated in m as 

tlit) = \ j ^dz [dlo{0) - 4o{e)] Tsu{t, z), (A.3) 


where 2 = cos9, are elements of Wigner’s d- 

matrix and 


Tsu{t,z) 


_4 ( 1 

3 \ml - s{t,z) 


1 


ml - u{t,z) 


(A.4) 


is the p-pole contribution in the crossed channels. The 
dimensionless coupling constants are defined as m 


9l — QpT^T^'t 92 — 9oJpTT 


(A.5) 


and the Mandelstam variables have the expressions 
s{t, z) = R{t) + K{t)z, u{t, z) = R{t) — K{t)z, (A.6) 




with 


TY}^ _(— *^777^ _ f 

^ , K{t) = 2q{t)p{t), (A.7) 


where q{t) and p{t) are defined below Eqs. ([2]) and (HI), 
respectively. 

The part t^{t) of the amplitude accounts for the rescat¬ 
tering in the direct channel. In the two-pion approxima¬ 
tion, unitarity gives 


disc [t^{t)V, ^(t)] =0, t> 4m^, 
where il{t) is the Omnes function 


= exp 


i r 

J4m^ t'it’ - t) 


(A.8) 


(A.9) 


is the energy-dependent p width defined in terms of the 
physical width Tp = 147.8 ± 0.9 MeV [s^, and 

777 h 

9{t) = (Ht) - Hml) -{t- ml)k\ml)) , 




T\/t 


2m„ 


(A.17) 


Choosing the subtraction point at to = ?T^pi the behav¬ 
ior of t^ (t) near t = m-p implies 

4 5152 
a = - ■ 


3 D{0)' 

Then, the reoresentation (IA.13I) is written finallv as 


which is analytic without zeros in the t-plane cut for t > 
4m^ and is normalized to 17(0) = 1. It can be written 
above the cut as 

+ ie) = t > 4 to^, (A.10) 

where 15^ (t) is the phase shift of the P wave of the elastic 
pion-pion amplitude. 

Since t\{t) is regular for t > 4mJ, from (IA.2I) and 
(IA.10|) it follows that 

disc [t)j(t)fl“^(t)] = —2it^(t)Im [^“^(t)]. (A.11) 


From the discontinuity one can reconstruct the function 
by _ means of a standard dispersion relation, written in 
as 


^{t) = a + 


t-to f°° fi(t')sin(5((t') 

71 JAm'^ - t) 

(A.12) 


in terms of the unknown subtraction constant a. Com¬ 
bined with (IA.2|) . this leads to 


t\t) 


n{t) 


' Ait) 

_n{t) 


H” d 


(A.13) 


i-h /■“ sin Slit') 

n J,^.jnit')\it'-to)it'-t) 


In , instead of the Omnes function il{t) a Gounaris- 
Sakurai parametrization [s^ was actually adopted, which 
is a reasonable approximation on the right hand cut 
where it is employed. Thus, 


nit) => GSit) = 


where 77(1) is written as 


rp(i) = 


i) - 

^ m ’ 

(A.14) 

1 - irripTpit). 

(A.15) 


(A.16) 


t\t) = 


m 


Ait)Dit) + g5i52 


(A.19) 


r 7 i(l')rp(l') 

Ti JatuI {f - ml)it' - t) 


dt 


As follows from (Q and (IA.4I1 . t\it) is real for t > 
4m^, which implies that the imaginary terms within the 
large parantheses compensate each other. Therefore, the 
phase of 1^(1) is equal to the phase of the Omnfe function 
1/17(1), i.e. to the phase shift ^((1). In this model, 1^(1) 
satisfies Watson theorem, and leads to a purely imaginary 
discontinuity © of the WTT form factor. As shown in 
[IIHIl, these properties are no longer valid in the more 
rigorous treatments of the amplitude. 

An obvious shortcoming of the model 0 is the fact 
that the amplitude 1/(1) was calculated in terms of a p- 
meson exchange neglecting the width of the p. In this 
approximation the w meson is actually stable since its 
mass is lower that the mass of pir pair. To improve the 
model, a straightforward procedure would be to include a 
finite width for the p poles in the denominators of (IA.4I) . 
We have adopted the prescription proposed in , where 
finite-width resonance exchange amplitudes with correct 
analyticity properties were obtained by replacing 

1 

"7-p - sit,z) 

and similarly for the u-channel contribution. The pole 
was replaced by a modified Breit-Wigner expression 
which automatically ensures the absence of singularities 
in the complex plane except for a right-hand cut. As 
suggested in [^, a reasonable choice for the spectral 
function cr(a:) is the imaginary part of the Breit-Wigner 
propagator: 


dx 


aix) 


JAm^ x-sit.zY 


(A.20) 


aix) 


_ rUpYpix) _ 

im?p - x)^ + m'^iTpix))^ ’ 


(A.21) 


with rp(a:) defined in (IA.16I) . In the limit of zero width, 
Tp —0, when cr(x) —?► irSix — rrip), the left side of (IA.20I1 
is recovered. 


In this relation 










































9 


By inserting the prescription (IA.20|) in (IA.3I) - (IA.4I) . the 
integration upon z = cos 9 can be performed exactly, 
leading to 


tlit) 


4gig2 

3 TT 



dx a{x) 


nz) 

K{t)’ 


where 


F{y) 


3 

4 


2g+(l 


g2)ln 


g + 1 ' 
y-1. 


and 


X — R{t) 

m 


(A.22) 


(A.23) 


(A.24) 


with R{t) and K{t) defined in (IA.7I) . 

The singularities of in the complex plane arise 

from the singularities of the function F{y) at y = ±.1 pro¬ 
duced by the logarithm, which depend parametrically on 
X. When x varies along the integration range in (IA.22I) . 
the singularities describes paths in the complex t-plane 
and in principle can overlap with the t-channel unitarity 
cut along the real semiaxis t > Arri^. The overlap can be 
avoided by a suitable prescription. In the present study 
we have adopted the prescription proposed in , which 
consists in adding to a small imaginary part, i.e. 

—>■ -I- le, with e > 0. With this prescription, we 

checked numerically that the singularities of {t) do not 
cross the unitarity cut in the t-plane. Moreover, the am¬ 
plitude t\{t) has no discontinuity across the line t > Am^, 
although it is no longer real on the unitarity cut. 

The points t = 4m^ and t = t±, i.e. the physical 
thresholds and the pseudo-threshold t_, require special 
attention since there the function K{t) defined in (IA.7I1 


vanishes. By using the asymptotic expansion 

F{y)r^- + -^ + ..., |y| > 1, (A.25) 

y 5g'^ 

and the decrease a{x) ^ l/x, we have checked explicitly 
that in the present model t\{t) is regular at these points. 

Since the amplitude t\{t) calculated from (IA.22I) has 
no discontinuity across the unitarity cut t > 4m^, the 
representation (jA.1911 remains valid. However, as men¬ 
tioned above, is complex for t > 4m^. Therefore, 

from (IA.19|) it follows that the phase of the amplitude 
t^{t) for t > 4m^ is no longer equal to the phase Sl{t) of 
the function 1/D{t). Watson theorem, which was valid in 
the original N/D model, is no longer valid now. More¬ 
over, the amplitude is not an analytic function of real 
type. These properties are satisfied of course by the ex¬ 
act solution fi{t) calculated in [iMl. 

It is useful to compare the simple improved N/D 
model presented here with the exact amplitude calcu¬ 
lated by solving numerically integral equations of the 
Khuri-Treiman type. An obvious feature of the N/D 
model is the lack of symmetry between the direct (t) and 
the crossed (s and u) channels. In fact, in the decay re¬ 
gion the dynamics in the three two-pion channels must 
be the same. In the Khuri-Treiman formalism, by iter¬ 
atively solving the relevant integral equation, the sym¬ 
metry between the three channels is gradually increased. 
This adjustment is not performed in the N/D approach, 
which has a rigid structure. However, by improving the 
description of the crossed channels in the frame of the 
N/D model, the main features of the exact partial wave 
amplitude fi(t), namely the failure of Watson theorem 
and the breakdown of the reality property, appear in a 
natural way. 
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